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Introduction: The role that Galactic cosmic rays 
(GCRs) and solar energetic particles (SEPs) play in the 
dynamic evolution of protosolar disks and the origin of 
our Solar System is a fundamental one. The GCRs are 
an important component of the interstellar medium 
(ISM), and even play a role in correcting the age de- 
terminations of some irons versus CAIs (calcium- 
aluminum inclusions) in meteoroids [1], Because CRs 
also are one of the energy transport mechanisms in a 
planetary nebula, the question of modelling their effect 
upon this broad subject is a serious topic for planetary 
science. The problem is addressed here. 

Relevance of the Modelling Problem: CRs are 
energetic nuclei that interact with Solar System matter, 
the physical nature and depth of collisional interaction 
being a function of energy. Chondrites, for example, 
are extraordinary mixtures of presolar and solar nebu- 
lar materials and asteroidal debris [1] that have under- 
gone numerous changes, some of which include nu- 
clear transmutation and heating induced by CRs. The 
original chemical and physical state of such material 
has thus been extensively obscured by complex proc- 
esses and isotopic change [2], some of which is re- 
flected in CR exposure histories [3]. The physics and 
chemistry are basically understood but the ambient CR 
environment is highly variable and usually unknown. 
This can be seen in the X-wind model of Shu et al. [4, 
5] and the bipolar jet model of Liffman & Brown [6], 

Modelling the effect of CRs during the formation of 
the Solar System is thus a daunting task. For example, 
one hypothesis is that 26 A1 may have been produced in 
the presolar molecular cloud due to irradiation by CRs 
[7, 2], Another is the complex nature of the proto- 
Sun’s journey through the local ISM (LISM), both past 
and present [8-22], Hence, CRs serve at least two roles 
during the origin of solar systems, both as a radionu- 
clide transmutation mechanism as well as for the trans- 
port of energy. These obscure the entire process of 
formation of what actually happened, yet is analyzed in 
present-day meteorites and extraterrestrial samples. 

Misconceptions: An additional problem regarding 
CR effects involves a serious misconception in the lit- 
erature. A number of authors [e.g. 23-24] have made 
claims that the astrophysical interaction of particles 
such as CRs in the ISM can be simulated using accel- 
erated-ion beams from Earth-based particle accelera- 
tors. It is true that the latter have been used for CR 
calorimeter calibrations [25] but this is not the same 
thing. The notion that accelerator beams can simulate a 
CR flux is false, and that fact has been demonstrated 


experimentally at CERN using CERF (CERN-EU 
High-Energy Reference Field) [26] embedded in the 
SPS (Super Proton Synchroton) beamline there. CERF 
attempted to moderate accelerator beams in an effort to 
simulate radiation fluxes and spectra at aircraft and 
space altitudes where neutrons and muons dominate. 
When CERF was adjusted to produce the neutron spec- 
tra, it failed for the muons - and vice versa. 

The reason is quite simple. GCRs and SEPs are 
single-event phenomena whereas accelerator beams 
have astronomical luminosities (comparatively speak- 
ing) and do not represent the same physics as CR inter- 
actions. CERF showed that the resultant muon fluxes 
during intranuclear cascades were wrong in comparison 
to CR-induced events. 

In summary, a comprehensive strategy for analyz- 
ing CRs in the protosolar disk is necessary and must 
involve satisfactory computational models using Monte 
Carlos, ISM transport codes, and dynamic modelling of 
the proto-Sun and heliosphere. 

CRs and Protoplanetary Disks: When the em- 
bedded young Sun appeared in a solar nebula forming 
from dense molecular clouds that made up the LISM of 
our Galaxy, a basic structure for the Solar System 
emerged. A solar wind appeared that created a bow- 
shock in the LISM as the proto-Sun rotated about our 
Galactic Center, and a heliosphere formed. The helio- 
sphere and its interaction with the LISM are the princi- 
pal regulatory mechanisms defining the ambient CR 
environment. Analytical studies are usually comprised 
of a hydrodynamic or magnetohydrodynamic (MHD) 
adaptation of two counterflowing fluids [18], the solar 
wind and the interstellar plasmas in a three-phase ISM 

[27] , Some energetic neutral particles are unaffected by 
the bowshock but later become photoionized CRs 
known as anomalous CRs (ACRs) [8, 10, 11, 15], 

The solar wind consists of energetic plasma with an 
embedded magnetic field that modulates the GCR flux. 
The modulation involved is usually based upon a time- 
dependent transport equation originating with Parker 

[28] that is derived from Fokker-Planck theory 
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where f(r, P, t) is the CR distribution function, P is 
rigidity, r is position, V is solar wind velocity, Q is 



some source function, and t is time. K represents the 
symmetric diffusion tensor, and one of the most rele- 
vant assumptions in CR modulation studies [8-22] is its 
component diffusion coefficients. 

The subject such as Eq. (1) has a long history in- 
volving MHD plasma theory (for what is known as the 
Navier-Stokes-MHD regime, X„«l) or the kinetic 
Boltzmann transport equation (known as the kinetic 
regime, K n » 1). Because MHD is a fluid theory (non- 
linear Navier-Stokes), it cannot be used for kinetic 
phenomena (Boltzmann equation). The difference is 
determined by the Knudsen number K„, with the transi- 
tion between being defined as K„ = 1 . 

Using these fluid and kinematic transport codes one 
can attempt to simulate the CR environment about the 
protosun for the X-wind [4, 5], the bipolar jet model 
[6], or any other as appropriate. At the present time, in 
contrast, our heliosphere is basically a low-density cav- 
ity carved out of the LISM by the solar wind [12], It is 
embedded in a warm, low-density interstellar cloud 
flowing through a local rest frame at ~18 km s' 1 . 

Galactic Propagation Models: Clearly the CR 
flux varies not only due to modulation by solar activity 
via the solar wind, but it also changes due to a variable 
heliosphere in response to its interaction with the dy- 
namic LISM. 

How can one address the basic question of CR- 
induced properties observed in cosmochronic archives 
such as ice cores or chondrites? The computational 
tools all exist, consisting of three basic categories: the 
plasma transport codes, the dynamic solar modulation 
codes, and the Monte Carlo interaction codes. The first 
two have been discussed above, while the third takes 
the modulated CR plasma flux as an input to induce 
changes in cosmogenic materials such as chondrites, 
primitive meteorites, interplanetary dust, asteroids, 
comets, and lunar rock or regolith. One of the best 
Monte Carlos is FLUKA used at CERN [29, 30], 

To the blend above must be added a newcomer that 
is seeing a number of important applications such as 
the Fermi Gamma-Ray Large-Area Telescope 
(GLAST). This is the existing Galactic Propagation 
(GALPROP) code available through Stanford Univer- 
sity and partially funded by NASA [31-34], The inter- 
pretation of current data from today’s CR and y-ray 
detecting satellites requires something beyond previous 
transport codes - which is one of GALPROP ’s primary 
goals. It is essential to upgrade many of the older exist- 
ing codes using measurements being made today in 
both astrophysics and CR physics [31, 35-37], 

Hence GALPROP represents an effort to incorpo- 
rate some of the latest spacecraft measurements such as 
EGRET [31] and other astrophysical data from BESS 
[35], ATIC [36], and Pamela [37] into the numerical 


propagation codes that deal with LISM and helio- 
spheric variations that influence CR flux. Similar modi- 
fications can prove useful in the study of CR contribu- 
tions to the protosolar disk and our own Solar System 
when the three basic code categories are adequately 
addressed. 
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